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The  Vaisala-Brunt  frequency  is  examined  for  data  obtained  at  an 
anchor  station  in  the  Florida  Current  for  a  ten  day  period.   There 
are  ninety-three  separate  S  T  D  stations  examined  to  yield  a  mean 
profile  of  the  Vaisala-Brunt  frequency  at  this  location. 

The  maximum  frequency  of  the  mean  profile  is  calculated  to  be 
0.0169  sec-1  and  occurs  at  a  mean  depth  of  100  meters.   This  cor- 
responds to  a  Vaisala  period  of  6.19  minutes. 

The  depth  at  which  the  maximum  frequency  occurs  varies  with 
time,  but  with  no  apparent  periodicity;  this  is  also  true  of  the 
magnitude  of  the  maximum  frequency.   This  randomness  may  be  the 
predominant  feature  of  the  Vaisala-Brunt  frequency. 

The  plot  of  frequency  versus  depth  for  each  station  reveals  fine 

structure  in  the  vertical  column,  showing  a  layered  structure.   The 

shallower  layers  appear  to  have  a  stronger  lamination  effect  on  the 
density  structure. 
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1,   INTRODUCTION 

Density  gradients  affect  the  motion  of  the  seas,   The  salinity, 
temperature  and  pressure  gradients  all  cause  the  sea  water  density  to 
vary  with  depth.   Most  of  the  density  gradient  is  caused  by  the  compres- 
sion of  the  water,  the  rest  by  the  variation  of  temperature  and  salinity 
with  depth  (Eckart,  1961).   Associated  with  the  density  gradients  are 
well  defined  density  boundaries  which  exhibit  certain  frequencies  and 
modes  of  oscillation.   Thus,  any  given  density  boundary  will  have  its 
own  normal  oscillating  frequency.   This  is  called  the  Vaisala-Brunt 
frequency  and  is  expressed  in  radians  per  second.   It  is  a  measure  of 
the  stability  of  oceanic  stratification  (LaFond  and  LaFond,  1966). 
The  distribution  of  this  frequency  is  one  of  the  most  important  dynam- 
ical features  of  the  ocean. 

The  concept  of  vertical  stability  in  the  ocean  was  expressed  by 

Sverdrup,  et  a_l. ,  in  1942  as  E  =  10"^   d o~  t  ,  and  appears  in  Richardson's 

dZ 
number  as   1   d  [p } .   Pollack  (1954)  defines  vertical  stability  as 

[/o  ]   dZ 
follows : 

1.  As  a  measure  of  the  restoring  force  (or  resistance  to  change) 
in  the  vertical  direction. 

2.  As  an  acceleration  per  unit  change  of  geometrical  depth;  the 
acceleration  may  be  expressed  per  unit  change  of  pressure 
also;  Hesselberg  used  this  approach. 

3.  As  an  expression  of  any  continuous  function  which  is  a 
parameter  of  the  resistance  to  deformation  of  the  fluid  column. 


4.    As  the  amount  of  energy  required  to  move  a  particle  of  unit 
mass  a  unit  displacement. 

In  definitions  2  and  3,  stability  is  a  function  of  the  difference 

between  an  observed  vertical  mass  distribution  and  that  of  a  neutral 

equilibrium  mass  distribution.   Therefore,  zero  corresponds  to  the 

neutral  position,  stability  is  indicated  by  the  positive  values  and 

instability  by  negative  values  (Pollak,  1954). 

According  to  Hesselberg,  stability  £  =  g  dP   ,  or  E  =  g  [d/3  -  (d/j  ], 

r    dZ  '        ?[dZ   (dZ)A] 

where  (d/y  is  the  density  gradient  of  neutral  equilibrium  or  the  mass 

(dZ)A 
distribution  in  a  column  of  homogeneous  water  (pressure=0,  salinity  and 

temperature=constants)  with  negligible  molecular  conductivity.   Bjerknes, 

et  al. ,  (1934)  noticed  the  adiabatic  density  gradient  (d/°  )  was  the  re- 

(dP  )A 

ciprocal  of  the  term  in  the  Lapalcian  expression  for  the  velocity  of 

2 
sound  C  =  (dP)  ,  (Pollak,  1954).   Therefore,  the  adiabatic  density 


(dA)A 

2 
gradient  can  be  replaced  by  1/C   (see  Appendix  A).   This  is  possible 

because  of  the  analogy  between  the  pressure,  p,  and  the  geometric  unit 

of  depth,  Z.   Pollak  equates  stability  in  the  geometric  units  as: 

E=  9 . 8 1  [d/  -  9.81/  ] ,  which  appears  similar  to  an  expression  devised 

P      [dP   "c2"] 
by  the  meteorologist  Vaisala  in  1928,  which  he  designated  N  ,  a  sta- 
bility expression  with  dimensions  of  frequency.   Eckart  mentioned  this 
expression  in  notes  of  a  seminar  given  in  1949.   According  to  Pollak, 
the  source  for  the  expression  was  a  lecture  by  Solberg  at  the  Univer- 
sity of  California  in  1946.   The  original  reference  could  not  be  located. 

We  can  examine  the  frequency,  N,  by  noticing  that  if  the  density 
of  a  fluid  is  changed  adiabat ical ly,  then  dP  =  C   d/  .  Vhere  dP  and  d/ 


are  the  changes  in  pressure  and  density,  and  C  is  the  speed  of  sound. 
If  dP  is  caused  by  changing  the  level  of  a  parcel  of  water  by  an 
amount  dZ  then  dP  =  -/*gdZ  if  we  assume  a  positive  Z  -  axis  up,  and  g 
is  the  acceleration  of  gravity. 

Elimination  of  the  term  dP  gives: 

1  dP  «  -g 

P     dZ    c? 

The  quantity  g/C2  is  identified  as  the  adiabatic  density  gradient  and 
is  that  part  of  the  actual  density  gradient  associated  with  compres- 
sibility effects.   Let  E  =  -  (1   d/°  +  £,  )  and  we  find  that,  generally 

0°   dZ    C2) 

speaking,  it  is  a  positive  value  in  the  seas,   The  quantity  N  in- 
troduced by  Vaisala  and  Brunt  is  equal  to  gE.   It  has  frequency 

dimensions  and  therefore  the  period,  T  =  2Y[  (Eckart,  1961). 

N 

In  considering  the  significance  of  the  coefficient  N,  Eckart 

(1960)  considers  the  oscillation  of  a  small  fluid  mass  displaced  from 

its  zero  position  and  allowed  free,  adiabatic  movement,  with  the 

positive  Z-axis  upward.   He  assumes  the  particle  of  fluid  is  in  a 

small  balloon-like  membrane  which  permits  changes  in  density  to  occur 

unimpeded,  and  thereby  maintain  the  outside  and  inside  pressures  the 

same.   The  displacement  causes  a  pressure  change  and  an  accompanying 

density  change  within  the  balloon,  so  that  d/°  inside  =  jdP.   The  out- 

C2 

side  density  change  will  be  : 

dP    outside  =  fe  dP,  where  r  is  the  displacement.   The  buoyant 
'  dZ        Y 

force  on  this  balloon  is  then  equal  to 

g  (d  P  outside  -  d  P    inside),     (1) 


which  is 


(df  -  PR)   =-/°N2£ ,  (2) 

(dZ     C?) 


where 

N2  =  -g  [1  d^  +  g  ]  (3) 

[7°   dZ    c2]  * 

If  the  frequency,  N,  is  real,  the  parcel  of  water  in  the  balloon  would 

oscillate  freely  at  the  frequency  N.   Instability  would  exist  if  N 

was  negative,  causing  the  ballon  to  leave  its  initial  position  (Eckart, 

1960). 

It  is  believed  that  N  is  the  limiting  frequency  for  the  free 

internal  waves  because  of  viscous  damping  (McLellan,  196.5).   Internal 

waves  can  exist  if  their  frequency  of  oscillation  is  greater  than  the 

inertial  frequency,  but  less  than  the  maximum  Vaisala  frequency, 

N   .   There  should,  therefore,  be  no  free  internal  waves  with  fre- 

quencies  above  N   ,  and  few  with  frequencies  above  the  local  Vaisala 
max 

frequency  in  a  particular  water  mass  (Cox,  1962).   In  examining  in- 
ternal waves  near  Bermuda,  Haurwitz,  ejt  al.,  (19.59)  used  the  Vaisala 
frequency.   They  associated  the  maximum  frequency  with  the  seasonal 
thermocline  and  the  second  maximum  found  at  a  deeper  depth,  with  the 
permanent  thermocline.   They  also  concluded  that  at  any  depth,  there 
should  be  no  frequency  larger  than  the  maximum  Vaisala  frequency. 
Hence,  this  may  be  designated  the  "cut-off"  frequency.   This  cut-off 
frequency  is  known  to  vary  with  time  and  space  and  must  be  computed 
for  each  vertical  density  distribution.   When  converted  to  period  the 
Vaisala  frequency  relates  to  the  natural  oscillation  governed  by  the 
stability  of  the  water  layers  (Lee,  1960). 

In  the  upper  four  thousand  meters  of  the  ocean  the  period  has  a 
wide  range,  from  one  minute  to  four  hours.  In  deeper  waters  the 


frequency  is  difficult  to  compute  because  the  terms  in  the  brackets 
of  the  equation  are  very  nearly  equal  and  of  opposite  sign.   At  the 
temperatures  and  pressures  found  in  the  deep  waters  the  values  of  the 
coefficient  of  thermal  expansion  are  doubtful,  and  at  best  are  extra- 
polated values.   These  extrapolated  values  create  large  uncertainties 
in  the  value  of  Vaisala  frequency.   Uncertainty  also  arises  from 
errors  in  measuring  small  salinity  gradients  (Eckart,  1961). 

An  investigation  of  the  Vaisala  frequency  in  the  Florida  Current 
was  conducted  using  pilot  data  collected  in  a  twenty-four  hour  period 
during  the  winter  of  1966  by  Dr.  Saul  Broida.   Two  locations  were  used 
with  stations  taken  alternately  between  them.   Density  parameters 
were  measured  in  situ  with  a  Bisset-Berman  Salinity-Temperature-Depth 
(STD)  system.   The  data  were  analyzed  at  fifteen  meter  depth  increments 
from  the  surface  to  seven  hundred  meters,   The  Vaisala  frequency  and 
period  were  computed  for  each  depth  increment,  and  plotted  versus 
depth. 

The  frequency  equation  as  stated  by  Eckart  was  used: 

n2  =-g  n      d/o  +  £  )   in  radians  per  second. 
(/>  dZ    C2) 

A  computer  program  was  devised  to  calculate  the  terms  of  the  equation 

and  solve  for  the  frequency  and  period,  where  T  =  2*TT  ,        Curve  fit- 

N 
ting  was  accomplished  by  a  cubic  spline  interpolation  (see  Appendix  B) . 

Analysis  of  the  resulting  pilot  data  disclosed  some  interesting 

features : 

1.   At  any  given  depth  the  frequency  varied  with  time.   Some 

periodicity  was  implied,  but  insufficient  data  points  did 

not  allow  any  definite  conclusions. 


2.  Minimum  period  values  occurred  at  the  boundary  of  the 
isothermal  layer  and  the  mixed  upper  layer. 

3.  The  shape  of  the  frequency  curves  indicated  a  multi-layered 
structure  of  laminations  of  density,  as  noticed  elsewhere 
in  the  oceans  by  Stommel  and  Fedorov  (1967),  and  others. 

4.  The  depth  of  maximum  frequency  varied  with  time. 


II.   METHODS 
A,    DATA  COLLECTION 

During  the  period  from  19  February  to  March  1,  1968  the  R/V 
P1LLSBURY  was  anchored  in  the  Florida  Current  in  460  fathoms  of  water 
at  Latitude  25-39.1  N  Longitude  079-48.0  W  (Fig.  1).   The  objective 
was  to  gather  temperature,  salinity,  and  depth  information  to  study 
density  fluctuations  in  the  Florida  Current,  and  analyze  the  Vaisala- 
Brunt  frequency  at  this  location.   This  cruise,  P6801,  was  the  longest 
continuous  anchor  station  completed  in  the  Florida  Straits. 

A  Bissett-Berman  Salinity-Temperature-Depth  probe  (STD)  system 
was  utilized  to  gather  the  JLn  situ  data.  Probe  casts  were  conducted 
at  two  hour  Intervals.   The  initial  lowerings  were  made  to  750  meters, 
but  an  unfortunate  wire  entanglement  later  limited  the  casts  to  a 
maximum  500  meter  depth.   The  STD  probe  utilizes  sensing  elements  which 
send  back  the  J^n  si^tu  data  as  a  frequency  signal  for  each  parameter. 
These  frequencies  are  recorded  on  graphic  plots  and  magnetic  tape. 
The  conversion  to  actual  values  of  depth,  temperature,  and  salinity 
Is  done  with  standard  conversion  tables  of  frequency  versus  the  para- 
meter measured.   To  maintain  the  calibration  and  accuracy  of  the 
system,  Niskln  bottle  samples  were  taken  for  comparison.   At  selected 
times  a  sample  was  taken  at  the  surface  and  one  at  the  maximum  depth 
along  with  the  probe  lowering.  The  salinity  and  temperature  readings  of 
both  instruments  were  compared  and  found  to  be  in  good  agreement. 


8 
Because  of  the  nature  of  the  Florida  Current  a  strong  surface 
current  prevailed  throughout  the  experiment,  approximately  4.2  knots. 
Technical  problems  and  rough  seas  limited  the  total  number  of  stations 
taken  to  ninety-three  during  the  ten  day  period. 

The  ship's  position  was  determined  by  visual  bearings  and  radar 
ranges  to  Fowey  Rocks  Light.  The  accuracy  of  the  anchor  position  is 
+  100  meters  in  an  East-West  direction;  +  400  meters  North-South. 

Some  gaps  do  exist  in  the  sequence  of  data  collection  which  are 
attributed  to  mechanical  failures  and  weather  conditions.   However, 
the  major  portion  of  the  stations  are  at  two  hour  intervals  with  the 
last  sequence  of  stations,  number  seventy-one  to  ninety- three,  being 
taken  at  hourly  intervals. 


FIGURE  1.    Chart  of  ship's  position 
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B.    DATA  ANALYSIS 

The  STD  data  from  the  cruise  were  in  analog  format  for  each 
station,  but  those  stations  taken  after  station  number  32  also  recorded 
the  parameters  on  magnetic  tape  via  a  digital  data  log  unit  (DDL). 
Those  stations  not  recorded  on  magnetic  tape  initially  were  digitized 
with  a  CALMA  converter.   This  is  a  manual  trace  of  the  temperature 
and  salinity  graphs,  which  records  the  horizontal  and  vertical  trace 
distances  on  magnetic  tape.   By  computer  conversion  the  distances  are 
changed  to  the  desired  parameters.   The  end  result  is  a  magnetic  tape 
record  of  all  the  station  data. 

The  STD  scanning  rate  was  set  at  0.1  seconds  giving  a  high  density 
of  points  on  the  traces  and  tapes.   To  select  the  most  representative 
value  of  a  parameter  at  a  given  depth  an  averaging  process  was  used. 
The  DDL  values  were  averaged  at  twelve  meter  intervals  from  the  sur- 
face to  the  maximum  depth  of  the  station  cast.   A  comparison  between 
the  averaged  values  and  the  actual  graphic  values  of  temperature  and 
salinity  at  given  depths  indicated  a  tolerance  of +  0.02c  G  for  tempe- 
rature and  +  0.03  o/oo  for  salinity.   The  averaging  process  provided 
the  best  values  for  analysis  for  the  total  cast. 

The  Temperature-Salinity  (T-S)  diagram  for  the  composite  of  the 
data  was  investigated  and  compared  to  the  historic  data  for  that  region. 
The  historic  data  was  taken  from  the  works  of  Parr  and  various  cruises 
conducted  by  the  Institute  of  Marine  Sciences.   Even  though  the  curves 
were  geometrically  coincident,  the  salinity  values  of  some  stations 
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appeared  to  be  lower  than  the  historic  data.   Although  the  spread  of 
the  T-S  curve  was  larger  than  might  be  expected,  the  relative  changes 
in  temperature  and  salinity  appear  adequate  for  analysis  of  the  gra- 
dients.  Thus,  these  relative  values  are  acceptable  for  the  investi- 
gation of  the  Vaisala  frequency  which  is  density  gradient  dependent. 

There  are  many  small  features  shown  on  the  temperature  traces 
which  are  probably  real,  but  the  salinity  traces  have  many  small  scale 
spikes  of  low  salinity  which  appear  false,  probably  caused  by  imperfect 
matching  of  the  time  constants  of  the  temperature  and  conductivity 
sensors  of  the  STD.   The  falsification  appears  as  fine  scale  variations 
in  salinity  at  the  points  of  strongest  vertical  temperature  gradient. 
The  deeper  the  traces  go,  the  more  the  salinity  trace  becomes  a  s^arp 
indicator  of  temperature  gradients. 

It  is  felt  that  many  of  these  spikes  are  compensated  for  in  the 
averaging  process  and  also  in  the  interpolation  of  points  while  curve 
fitting. 
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C.    COMPUTATION 

Once  the  station  data  was  compiled  in  terms  of  temperature, 
salinity,  sigma  -  t,  and  specific  volume  anomaly,  for  each  twelve  meter 
increment  of  depth,  the  parameters  for  solving  the  Vaisala-Brunt  fre- 
quency equation  were  avialable.   The  value  for  the  acceleration  of 
gravity  was  assumed  constant  at  9.7893  meters  per  second  at  Latitude 
25.7  N.   Sound  speed  was  computed  at  each  depth  increment  by  the 
following  equation: 

C  =   141,000  +  421  t  -  3.7  t2  +  110  S  +  0.018  d,  in  (cm/sec),  (4) 
where 

t  =   temperature  (°C),  S  =  salinity  (o/oo),  d  =  depth  (cm). 
These  values  of  speed  were  converted  to  meters  per  second. 

The  _in  situ  density  was  computed  and  used  in  determining  the 

quantity  d/3  .   The  equation  used  to  find  the  density  was: 
dZ 

OC  (35,0,P)+S  =  OC  in  situ    (cnvVqm),  (5) 

where 

OC(35,0,P)=  standard  specific  volume  at  depth,  and  S  =  specific 
volume  anomaly.   The  reciprocal  of  the  specific  volume  _in  situ  is  the 
in  situ  density,  P  .   The  standard  specific  volume  was  computed  by 
linear  interpolation.   The  values  decreased  linearly  with  depth  by 
45  X  10" 7  cm3/  gm/  meter. 

The  values  of  density  jLn  situ  were  plotted  versus  depth  for  each 
station  and  a  cubic  spline  interpolation  fit  was  made  to  the  points, 
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(Fig.  2).   Using  interpolated  values  from  these  curves  the  slopes, 

dJ2_,   were  found  at  five  meter  intervals,  as  were  the  values  of  density 

dZ 

at  those  depths. 

The  frequency  equation  was  solved  using  these  interpolated  values 

in  the  equation: 

N2  =  -g  (1  dP    +  a  ).  (6) 

(r    dZ    C2) 

Frequency  was  converted  to  period  by  T  =  2*TT   minutes.   The  frequency 

N(60) 
was  plotted  versus  depth  and  analyzed,  (Fig.  3).   A  mean  frequency 

curve  was  also  plotted  for  the  entire  cruise  data  by  averaging  the 

interpolated  values  of  frequency  at  five  meter  increments  for  the 

ninety-three  stations  (Fig.  4).   Where  negative  values  of  N  were 

encountered,  they  were  eliminated  from  the  average  to  present  only 

positive  values  of  frequency.   These  negative  values  are  interesting 

in  that  if  they  really  exist  they  could  represent  instability  in  the 

laminar  structure  of  the  water  column,  but  they  may  be  caused  by 

erroneous  data  readings  as  a  result  of  time  lags  in  the  sensors,   It 

would  be  interesting  to  investigate  these  in  future  studies.   The 

historic  data  from  the  Florida  Current  indicate  strong  stability  in 

the  area  of  the  experiment;  therefore,  it  is  felt  that  these  negative 

values  might  really  be  neutral  stability  points,  approaching  N=0. 

Because  they  were  few  and  variable  in  depth,  it  is  felt  that  the 

averaging  process  was  not  prejudiced  by  their  elimination,  nor  was  the 

investigation  of  the  frequency  in  depth  and  time.   When  possible,  depths 

were  chosen  where  a  minimum  or  no  negative  values  of  N  were  noted 

during  analysis. 

To  help  diminish  errors  in  using  Vaisala  frequency  as  a  measure 

of  stability,  small  depth  increments  were  chosen  for  analyzing  gradients 
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FIGURE  2.    Density  versus  depth  profiles 

(6  stations) 
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FIGURE  3.    Frequency  versus  depth  profiles 

(6  stations) 
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FIGURE  4.    Mean  frequency  versus  depth  profile 
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in  density.   Hence,  the  five  meter  increment  for  interpolated  values 
of  density.   According  to  Roden  (1968)  it  is  believed  that  any  fun- 
damental oscillations  would  differ  from  the  Vaisala  frequency  by  about 
twenty-five  per  cent  for  this  vertical  displacement.   In  studying  the 
Vaisala  frequency  in  the  Gulf  of  Alaska,  Roden  found  that  the  non- 
linear cases  showed  stability  oscillations  no  longer  symmetrical  about 
the  mean,  as  is  the  Vaisala  oscillation,  and  the  fundamental  fre- 
quency of  oscillation  differed  from  the  Vaisala  frequency. 

In  computing  the  standard  deviation  and  variance  of  the  fre- 
quency for  the  entire  period  of  observations,  100,  245,  and  350  meters 

were  selected  as  representative  depths.   They  were  so  chosen  because 

9 
they  contain  a  minimum  of  negative  N  values,  which  were  to  be  ex- 
cluded; they  included  the  mean  depth  of  the  maximum  frequency  as 
compiled  in  the  mean  frequency  curve;  and  they  are  spaced  far  enough 
apart  to  give  an  in  depth  view  of  the  frequency  on  a  statistical 
basis . 

From  each  station  the  maximum  frequency  was  gathered  and  plotted 
versus  time  (Fig.  5),  and  so  was  the  depth  at  which  the  maximum  fre- 
quency occurred  (Fig.  6).   These  were  visually  inspected  for  some 
recognizable  periodicity,  but  none  was  readily  apparent.   In  fact, 
the  random  nature  of  the  function  is  the  primary  observation, 

From  the  mean  profile  of  Vaisala  frequency  a  maximum  value  of 
0.0169  sec"   was  obtained,  which  corresponds  to  a  period  of  6,19 
minutes.   The  depth  of  the  maximum  frequency  was  100  meters.   The 
computed  standard  deviation  of  depth  for  maximum  frequency  was  +  15,2 
meters,  with  a  relative  dispersion  (coefficient  of  variance)  of  16.6 
per  cent. 
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The  standard  deviation  of  the  frequency  at  100  meters  was  + 
0.0037  sec"  ;  with  a  relative  dispersion  of  22.8  per  cent.   At  245 
meters  the  standard  deviation  was  +  0,0030  sec"   and  the  relative 
dispersion  was  40-5  per  cent.   For  350  meters  the  standard  deviation 
was  +  0.0031  sec"   with  a  relative  dispersion  of  60,9  per  cent. 


23 


FIGURE  5.    Maximum  frequency  versus  time 
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FIGURE  6„    Depth  of  maximum  frequency 

versus  time 
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III.    RESULTS 
A.    DISCUSSION 

Results  of  the  experiment  discussed  in  this  thesis  lead  to  the 
observation  that  the  Vaisala-Brunt  frequency  in  the  Florida  Current 
near  the  axis  of  the  current  is  a  complex  feature  of  stability  with 
no  apparent  periodicity  with  respect  to  depth  or  time-   An  examination 
of  each  plot  of  frequency  versus  depth  shows  that  the  maximum  fre- 
quency occurs  in  the  thermocline  where  the  mean  density  gradient  is 
greatest  and  decreases  above  and  below  this  level  as  the  water  becomes 
more  homogeneous.   The  smallest  frequency  values  are  observed  in  the 
deeper  layers  and  near  the  surface. 

The  weather  during  the  data  collection  phase  was  particularly 
stormy  at  the  anchor  station.   This  may  account  for  the  well  mixed 
upper  layer  and  fairly  abrupt  thermoclines  recorded.   Thus,  we  ex- 
perienced a  three  layer  oceanic  structure,  with  relatively  homo- 
geneous water  above  and  below  the  thermocline  layer,   The  plots  of 
frequency  versus  depth  show  a  sharp  maximum  frequency  at  or  just  below 
the  thermocline  boundary  with  smaller  values  of  N  elsewhere.   The 
thermocline  thickness  appears  to  be  between  100-150  meters. 

It  is  observed  that  the  STD  data  provide  continuous  profiles  of 
the  thermohaline,  and  therefore  density,  structure.   The  instrument 
makes  obsolete  the  older  methods  of  interpolating  between  widely  spaced 
sample  points.   The  STD  reports  discontinuities  varying  in  magnitude, 
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time,  and  space  as  the  rule  rather  than  the  exception.   In  the  plots 
of  frequency  versus  depth  (Fig.  3),  series  of  layers  are  observed 
clearly  with  sharp  boundaries  separating  each  from  its  neighboring 
layers,  and  these  exist  throughout  the  water  column.   The  gradients 
or  boundaries  between  the  layers  are  often  thin  and  sharp.   The 
shallower  layers  show  stronger  lamination  effects  in  the  density 
structure  than  do  the  deeper  depths.   These  are  obvious  observations 
of  the  data  collected  and  are  only  the  beginning  of  exploration  into 
the  recognition  of  this  fine  structure.   Stommel  and  Fedorov  (1967) 
present  some  interesting  ideas  on  the  dynamic  processes  involved  in 
this  fine  structure  and  possible  gradient  effects  on  vertical  mixing 
and  mean  turbulent  transport.   Data  similar  to  that  obtained  for  this 
thesis  could  be  analyzed  in  the  light  of  their  studies  and  add  much 
to  the  study  of  the  Florida  Current. 

It  is  recognized  that  the  Vaisala-Brunt  frequency  plays  a  major 
role  in  the  study  of  internal  waves,  therefore  its  distribution  with 
depth  in  the  Florida  Current  is  of  interest,   As  noted,  it  is  a  very 
complicated  function  of  depth.   It  appears  that  the  complex  structure 
occurs  because  of  temperature  and  salinity  perturbations. 

The  average  plot  of  frequency  versus  depth  presents  a  representa- 
tive trace  of  the  Vaisala-Brunt  frequency  for  the  Florida  Current 
during  the  period  of  the  cruise.   The  small  deviations  indicate  that 
the  average  frequency  figure  might  be  useful  in  future  calculations. 

B.    RECOMMENDATIONS 

In  looking  ahead  to  future  work,  it  is  recommended  that  this  type 
cruise  be  conducted  again  in  the  summer  months,  preferably  to  obtain 
a  seasonal  variation  of  the  density  structure  and  see  the  effect  on 
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the  Vaisala-Brunt  frequency.   More  emphasis  should  be  placed  into 
examining  the  fine  structure,  however,  because  this  may  lead  to  a 
better  understanding  of  how  the  density  structure  varies,  and  why. 
It  may  also  offer  an  understanding  of  the  Vaisala  frequency  randomness 
of  variation. 

Project  MIMI,  a  study  of  acoustic  signal  propagation  across  the 
Florida  Straits,  collected  acoustic  data  simulaneously  with  our  col- 
lection of  density  data.   An  examination  of  these  data  might  provide 
an  insight  into  the  internal  wave  configurations  of  the  Florida  Current, 
or  produce  a  means  of  measuring  Vaisala-Brunt  frequency  directly,  if 
a  definite  relationship  exists  between  the  variations  of  N  and  the 
acoustic  variations  received  in  Bimini.   It  offers  an  excellent  op- 
portunity to  compare  physical  parameters  with  acoustic  data  taken 
simultaneously  over  an  extended  period,  along  the  acoustic  axis. 

Any  future  investigations  of  the  Vaisala-Brunt  frequency  in  the 
Florida  Current  and  its  distribution  will  enhance  the  dynamical 
investigations  of  this  region. 
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IV.   SUMMARY 

The  Vaisala-Brunt  frequency  was  calculated  and  investigated  for 
a  given  location  in  the  Florida  Current  near  the  axis  of  the  current. 
The  analysis  of  data  from  ninety-three  separate  STD  stations  gives  a 
composite  profile  of  the  Vaisala-Brunt  frequency  versus  depth  for  this 
region.   The  experiment  was  designed  to  obtain  density  information  for 
an  investigation  of  the  fluctuations  of  the  frequency  as  a  function 
of  time  and  depth.   The  results  of  the  investigation  are  station 
profiles  of  the  frequency  versus  depth,  a  mean  profile  for  the  aggregate 
of  stations,  a  plot  of  maximum  frequency  versus  time  and  the  depth  of 
N  max  versus  time.   The  standard  deviation  and  variance  were  computed 
at  three  selected  depths  for  frequency,  and  for  the  depth  of  N  max. 

The  investigation  reveals  that  the  frequency  is  a  very  complex 
function  of  depth  and  time  which  appears  to  vary  in  time,  space,  and 
magnitude  in  a  random  manner;  typical  of  many  oceanic  processes.   This 
randomness  of  the  function  is  the  most  significant  observation.   A 
representative  figure  for  N  max  is  presented  along  with  the  depth  at 
which  it  occurs  for  possible  use  in  future  work  related  to  internal 
waves,  acoustics,  or  density  stability  in  the  Florida  Current, 

CONCLUSIONS 

The  mean  value  of  maximum  Vaisala-Brunt  frequency  at  this  location 
in  the  Florida  Current  is  0,0169  sec"   with  a  corresponding  period  of 
6. 19  minutes . 
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The  average  depth  of  occurrence  of  the  maximum  Vaisala-Brunt 
frequency  is  LOO  meters  with  a  standard  deviation  of  +  15,2  meters 
and  a  coefficient  of  variance  of  16,6  per  cent,  or  relative  dispersion. 

The  standard  deviation  of  the  mean  Vaisala-Brunt  frequency  at 
100  meters  is  +  0.0037  sec"1-  with  a  coefficient  of  variance  of  22-8 
per  cent.   At  245  meters  the  standard  deviation  decreases  to  +  0.0030 

sec"   and  the  coefficient  of  variance  increases  to  40.5  per  cent.   At 

_  i 
350  meters  the  standard  deviation  is  +  0.0031  sec   with  a  coefficient 

of  variance  of  60.9  per  cent.   All  the  coefficients  of  variance  are 

increasing  with  depth,  with  the  standard  deviation  fairly  constant  with 

depth.   The  deviation  is  slightly  higher  within  the  thermocline;,  near 

maximum  frequency. 

The  average  plot  of  frequency  versus  depth  is  shown  in  figure  4 
and  is  representative  of  the  Vaisala  frequency  in  depth  for  this 
location. 

The  depth  at  which  the  maximum  frequency  occurs  varies  with  time 
as  does  the  magnitude  of  maximum  frequency,  however,  no  apparent 
periodicity  is  noticed.   The  randomness  of  this  feature  may  be  its 
most  important  characteristic. 
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APPENDIX  A.   VELOCITY  OF  SOUND  TERM  IN  THE 
VAISALA- BRUNT  EQUATION 
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From  Hook's   Law:  STRESS   =  K     X     STRAIN. 


E  =  Bulk  Modulus,    in  dimensions   of  pressure   FORCE  or  MASS  x  ACCELERATION 

AREA  AREA 

P  ■  density,  in  dimensions  of  MASS 

VOLUME 

so  that   E  =  ML    =   M    and  P    =  M  ,  (1) 

T^L2     LT2"  \? 

E  =  M   x  iP   =  L2   =   (  L  )2  (2) 

/o   ltZ     m      ^2"     (  t  ) 

which  is  in  velocity  units  squared. 


According  to  Newton,  the  speed  of  propagation  of  waves 

C2  =  E  ,  (3) 

7° 

substituting  pressure  P,  for  E,  we  obtain  the  same  results,  so  that 

P   =  C2,  (4) 

7 

and  for  an  adiabatic  change 


(dP  )   =  C2.  (5) 

(d/°)A 
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APPENDIX  B.    CUBIC  SPLINE  INTERPOLATION 
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A  new  approach  based  on  an  old  idea  is  used  to  interpolate  a  set 
of  data.   It  is  based  on  the  idea  of  using  a  different  expression  in 
each  sub- interval  [X.,  X^+jJ  an^  auxiliary  conditions  of  smoothness 
at  the  joints  instead  of  a  simple  closed  form  expression,  valid  over 
the  entire  interval.   It  was  developed  because  the  other  methods  (e.g., 
least  squares)  failed  to  give  good  results  under  certain  conditions; 
at  the  two  extremes  of  a  set  of  data  or  could  not  provide  a  good 
evaluation  of  the  first  derivative  along  the  computed  polynomial  fit. 

The  basis  for  all  interpolation  schemes  is  that  given  a  discrete 
set  of  n  +  1  points  an  nth  degree  polynomial  can  be  passed.   An  inter- 
polating polynomial  of  degree  as  high  as  the  number  of  data  points 
would  allow  can  be  obtained.   But  this  can  be  impractical  sometimes. 
Some  of  the  problems  encountered  with  interpolating  schemes  like  the 
least  square  method  are: 

1.  A  set  of  data  having  points  which  are  far  apart  from 
each  other,  will  fail  to  give  a  good  approximating 
polynomial,  since  the  conditions  for  smoothness  are 
very  fragile,  in  most  cases. 

2.  If  the  change  in  the  slopes  from  point  to  point  are 
very  sharp,  the  resultant  polynomial  produces  gross 
errors,  since  this  rate  of  change  of  slopes  has  a 
direct  effect,  not  only  where  the  change  occurs  but 
also  along  the  entire  interval. 

3.  Theory  states  that  the  higher  the  degree  of  computed 
polynomial  the  better  the  approximation.   Again, 
experience  shows  that  one  could  go  only  so  high  on  the 
degree  of  the  interpolated  polynomial  before  the  errors 
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pile  up  making  the  whole  thing  frustrating,  when  a 
considerable  amount  of  data  points  are  available* 
Experience  shows  that  polynomials  higher  than  13th 
degree  are  impractical.   But  this  depends  on  the 
distribution  of  the  data  points. 
All  this  led  to  the  development  of  the  Spline  interpolation  which 

is  based  on  the  idea  of  passing  a  low  degree  polynomial  through  a  few 

points  at  a  time  making  use  of  conditions  of  continuity  and  existence 

of  first  and  second  derivative. 

A  cubic  polynomial  is  interpolated  in  each  sub  interval  [X|"l,^\], 

i  =  l,...n.   It  is  accomplished  with  the  following  conditions  met, 

1.  Continuity  of  the  function 

2.  Continuity  of  first  derivative 

3.  Continuity  of  second  derivative. 

Due  to  the  nature  of  the  computations,  data  which  is  equally  spaced 
provides  short  cuts  to  the  procedure.   With  this  in  mind  the  data  for 
this  experiment  was  chosen  at  equal  intervals  of  depth. 

For  a  mathematical  treatment  of  the  Cubic  Spline  Interpolation 
scheme  the  reader  is  referred  to  "Approximation  of  Functions"  by 
Carl  de  Boor  (1965).   Elsevier  Publishing  Co.,  New  York. 
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APPENDIX  C,    STATISTICAL  FORMULAS 


n     =     number  of   points 
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MEAN  FREQUENCY 


X  =  X_ 

n 


(1) 


STANDARD  DEVIATION 


\ 


(X)' 


(2) 


VARIANCE 


V  =   S 


CO 


COEFFICIENT  OF  VARIANCE  (relative  dispersion) 

C  =   S   =   standard  deviation 


X 


mean 


(4) 
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